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■ Abstract. The 2:1 orbital resonances of the GJ 876 system can be eas- 
ily established by the differential planet migration due to planet-nebula 
interaction. Significant eccentricity damping is required to produce the 
observed orbital eccentricities. The geometry of the GJ 876 resonance 
configuration differs from that of the lo-Europa pair, and this difference 
is due to the magnitudes of the eccentricities involved. We show that a 
large variation in the configuration of 2:1 and 3:1 resonances and, in par- 
ticular, asymmetric librations can be expected among future discoveries. 

> 

^ - 

Q"^ ■ 1. Introduction 

o : 

CN . Marcy et al. (2001) have discovered two planets in 2:1 orbital resonances about 

the star GJ 876. This and possibly two other resonant pairs, one in 2:1 resonances 
. about HD 82943 (Mayor et al. 2001; Gozdziewski & Maciejewski 2001) and the 

other in 3:1 resonances about 55 Cnc (Marcy et al. 2002), and the ease of capture 
Q ■ into these resonances from nebula induced differential migration of the orbits 

j-j \ (e.g., Lee &: Peale 2002), mean that such resonances are likely to be ubiquitous 

■ among extrasolar planetary systems. 
Cd ■ 

> ' 

V-j ■ 2. The GJ 876 System 

■ A dynamical fit to the radial velocity data (e.g., Laughlin & Chambers 2001) 
places the GJ876 planets on coplanar orbits deep in three resonances at the 2:1 
mean-motion commensurability. The libration of the mean-motion resonance 
variables, = Ai — 2A2 + wi and 6*2 = Ai — 2A2 + and the secular resonance 
variable, 9^ = w\ — wi = Q\ — 62, about 0° (where Ai^2 are the mean longitudes 
of the inner and outer planet and wi^2 are the longitudes of periapse) differs 
from the familiar geometry of the 2:1 resonances between Jupiter's satellites lo 
and Europa, where 9i librates about 0° but O2 and ^3 librate about 180°. We 
have found that this difference is mainly due to the magnitudes of the orbital 
eccentricities involved ((ei) < 0.005 for the lo-Europa system, but (ei) 0.26 
for GJ 876). Stable simultaneous librations of the resonance variables require the 
equality of the average periapse motions, {wi) = {^2). When the Cj are small, 
the Wj are dominated by the lowest order resonant terms in the disturbing 
potential, and (tbi) = {W2) requires 0i « 0° and 62 ~ 180°. On the other 
hand, when the ej are large, the dominance of higher order terms means that 
{wi) = {W2) requires both 61 and 62 ~ 0° (see Lee & Peale 2002 for details) and 
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Figure 1. Evolution of the GJ876 planets into 2:1 resonances for a 
calculation where the outer planet is forced to migrate inward and there 
is no eccentricity damping (see text for the curves labeled K = 100). 
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Figure 2. Evolution into 2:1 resonances for a calculation with 
mi/m2 = 1.5 and (mi +m2)/mo = 10"^. 



for some systems with masses different from those in GJ 876, 9i and 62 far from 
either 0° or 180° (see § 3). 

The likely origin of the resonances in the GJ 876 system is the differential 
planet migration due to planet-nebula interaction. The planets are sufficiently 
massive to open gaps in the nebula disk surrounding the young GJ 876 and to 
clear the disk material between them (Bryden et al. 2000; Kley 2000), and the 
resulting planet-nebula interaction typically forces the outer planet to migrate 
inward on the disk viscous timescale (e.g.. Ward 1997). Figure 1 shows the 
evolution of the semimajor axes aj, eccentricities ej, and resonance variables 9i 
and 62 for a calculation where the outer planet of the G J 876 system is forced 
to migrate inward and there is no eccentricity damping. Both 9i and 62 are 
captured into libration about 0°, but the ej quickly exceed the observed values 
(dotted lines) with continued migration within the resonance. Therefore, unless 
resonance capture occurred just before nebula dispersal, eccentricity damping is 
required to produce the observed ej of the G J 876 system. Eccentricity damping 
of the form 62/62 = —Xl 02/02] causes the eccentricities to reach constant values 
that decrease with increasing K. Figure 1 also shows the evolution of the ej 
for a calculation with K = 100, which is chosen to yield the observed ej. In 
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Figure 1 we have neglected the outward migration of the inner planet due to 
the torques induced by any nebula left on the inside of the inner planet, as this 
docs not qualitatively change the evolution except to require less eccentricity 
damping {K fs 10) to arrive at the observed eccentricities (Lee & Peale 2002). 



3. Diversity of 2:1 Resonance Configurations 

As the example in Figure 1 shows, differential migration of planets through 
capture into a set of resonances and subsequent evolution without eccentricity 
damping cause the system to pass through a continuous sequence of stable li- 
bration states, whose configurations depend on the values of the eccentricities 
ej and the mass ratios mi/m2 and (mi + m2)/mo (where mo is the mass of 
the primary). Where a system with a given set of parameters is left along this 
sequence depends on the termination of the migration through nebula dispersal, 
or much more probably, the magnitude of the eccentricity damping. 

For (mi+m2)/mo ~ 10~^ (a factor of 10 smaller than in the GJ876 system), 
mi/ 1712 ^ 1, and very small initial ej, a system is captured initially into a lo- 
Europa like configuration, with 6i librating about 0° and 62 librating about 
180°. Continued migration forces ei to larger values and 62 from increasing 
to decreasing until the system passes smoothly over to the GJ 876 configuration 
with all 9j librating about 0°. Then both ei and 62 increase smoothly as in Figure 
1 until instability ensues at very large ej. For 1 ^ mi/m2 ^ 2.5 (see Fig. 2), 
the passage from the lo-Europa configuration to the GJ 876 configuration is as 
above, but when the increasing ei reaches ~ 0.1-0.3 (depending on mi/m2), the 
libration centers depart by tens of degrees from 0°, leading to stable libration 
of the 0j far from either 0° or 180°. This asymmetric libration can persist 
for a wide range of ei, but the system returns to the configuration with all 9j 
librating about 0° for sufficiently large ei. If mi/ 1712 ^ 2.5 (see Fig. 3), the 
resonance variables never reach the GJ876 configuration after leaving the lo- 
Europa configuration. Within each type of evolution, the details of the path 
depend only on mi/m2 and not on the total planetary mass, provided mj are 
not too large. If (mi + m2)/mo ~ 10~^ as in the GJ876 system, fairly large 
Cj are generated before resonance encounter and the lo-Europa configuration 
expected at small ej is skipped (see Fig. 1). 
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Figure 4. Evolution into 3:1 resonances for a calculation with the 
same masses as in Fig. 3. 



4. The 55 Cnc System and 3:1 Resonance Configurations 

For the inner two planets of the 55 Cnc system [with minimum (mi +7712) /mo ~ 
10~^] and the expected range of migration rates, capture into the 3:1 resonances 
is sensitive to the initial conditions and migration rate. Figures 3 and 4 show the 
results of two calculations that are different only in the initial orientation of the 
orbits, and the former is captured into 2:1 resonances while the latter is captured 
into 3:1. As shown in Figure 4, there arc also 3:1 resonance configurations with 
asymmetric librations. The "normal" behavior for small eccentricities where 
the 3:1 resonance variables (see Fig. 4 for definitions) librate about 0° or 180° 
changes to librations about angles far removed from cither of these as the system 
is forced deeper into the resonances with continued migration. 
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